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Non-Markovian decoherence in the adiabatic quantum search algorithm 
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We consider an adiabatic quantum algorithm (Grover's search routine) weakly coupled to a rather 
general environment, i.e., without using the Markov approximation. Markovian errors generally 
require high-energy excitations (of the reservoir) and tend to destroy the scalability of the adiabatic 
quantum algorithm. We find that, under appropriate conditions (such as low temperatures), the 
low-energy (i.e., non-Markovian) modes of the bath are most important. Hence the scalability of the 
adiabatic quantum algorithm depends on the infra-red behavior of the environment: a reasonably 
small coupling to the three-dimensional electromagnetic field does not destroy the scaling behavior, 
whereas phonons or localized degrees of freedom can be problematic. 

PACS numbers: 03.67.Pp, 03.67.Lx, 03.67.-a, 03.65.Yz. 



In contrast to (conventional) sequential quantum al- 
gorithms in which the solution to a given problem is 
obtained by measuring a final state after a sequence of 
operations (gates) acting on an initial state, the underly- 
ing idea of adiabatic quantum computation is to encode 
that solution in the ground state of a suitably designed 
Hamiltonian instead []|. Since the usual relaxation into 
that ground state typically takes an exponentially long 
time for a complicated and strongly coupled Hamilto- 
nian (as the system is very likely to be trapped in a local 
minimum), the desired ground state is reached by ex- 
ploiting the adiabatic theorem: Given a time-dependent 
Hamiltonian H{t) with its instantaneous energy eigenba- 
sis H{t)\Ek{t)) = Ek{t)\Ek{t)) , the actual quantum state 
of the system stays near the ground state |-Eo(0) 
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provided that the evolution is slow enough, i.e., that the 
first-order corrections of this adiabatic expansion (second 
term on the r.h.s.) are small. Note that the dynamical 
phase i^fc(t) — i /q drEkir) and the (geometrical) Berry 

phase 7fc(0 = i dT{Ek{T)\Ek{T)) have been absorbed 
into \Ekit)) for brevity. Therefore, by slowly evolving the 
quantum system from an initial Hamiltonian with an easy 
to prepare ground state (the initial state) into the final 
problem Hamiltonian, the system ends up in (or close 
to) the desired ground state (encoding the sought-after 
solution). For infinite running times, the system would 
perfectly stay in its ground state as the non-adiabatic 
corrections become arbitrarily small. But finite running 
times still allow for a reasonably high fidelity of the final 
ground state - the minimum runtime T is then a measure 
for the computational complexity of the algorithm. 

One of the major motivations for adiabatic quantum 
computation (and the main advantage in comparison to 
sequential quantum algorithms) is the relative robustness 
of the ground state against decoherence caused by the 
inevitable coupling to the environment - which makes it 



a very promising candidate for an experimental realiza- 
tion. However, so far there are only a few quantitative 
investigations of the impact of decoherence on adiabatic 
quantum computers. In addition, these investigations 
predominantly rely on master equations 0,01 or superop- 
erators (which are both based on the Markov approx- 
imation) or classical (Markovian) noise Q . But such 
processes require short time-scales and, in turn, large en- 
ergy scales - which is a cause for concern that the Markov 
approximation might not be suitable for adiabatic quan- 
tum computation, in which the system and the environ- 
ment are supposed to have a rather low temperature. In 
order to bridge this gap, we consider an adiabatic quan- 
tum algorithm (Grover's search routine) weakly coupled 
to a rather general environment and calculate the im- 
pact of decoherence directlv, i.e., without resorting to 
the Markov approximation 

The quantum algorithm we are going to study is an 
adiabatic version of Grover's search routine Q defined 
by the Hamiltonian (acting on an n-qubit system) 

H{t) = [l-s{t)\{t-\i^^){M)+<t){t~\w){w\), (2) 

with the initial ground state I^/jq) = Y^^Zq \x) /^fN being 
the coherent superposition of all numbers x and the fi- 
nal ground state \w) € {|a;)}^Q^ denoting the problem's 
solution. It solves the problem of finding a marked item 
w in an unstructured list of length iV = 2" and was first 
introduced in the framework of sequential quantum com- 
putation Note that not the marked state \w) itself 
but only the projector \w) {w\ onto this state is used in the 
Hamiltonian (similar to a quantum oracle: not knowing 
the solution, but being able to recognize it). We choose 
the adiabatic quantum search algorithm because its pa- 
rameters such as the instantaneous energy eigenbasis and 
the fundamental energy gap can be calculated analyti- 
cally - but we expect our results to reflect the typical 
behavior of general avoided level crossings. 

The time-dependence is governed by a strictly mono- 
tonic function s{t) that interpolates between the initial 
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and the final Hamiltonian with s(0) = and s{T) = 1, so 
that T is the algorithm's running time. The essential dy- 
namics take place in a two dimensional subspace spanned 
by {{w), {ipo}}- After a Gram-Schmidt orthogonalization, 
the basis we use for our calculations reads {\w), Iw-^)} 
with \w^) = (//ViV'o) - \w))/y/N -1. The fundamen- 
tal (non-vanishing) energy gap between the ground state 
and the first excited state 

AE = ^l + is{t)[l-sit)] (^-l) (3) 

exhibits a minimum at s = 1/2 of Ai?min — 1/%/^- 
The total running time necessary for keeping the non- 
adiabatic corrections small depends on the interpolation 
dynamics s{t). A constant velocity s = const yields 
the same scaling T = 0{N) as the classical (brute- 
force) search, cf. Eq. QJ. Adapted interpolations such 
as s oc (AE')^ or s oc Ai?, on the other hand, recover the 
quadratic speed-up T = 0{^/N) or T = 0(//V IniV), 
respectively 0. We shall use these three specific cases 
of s(t) for investigating the impact of a coupling to an 
environment. 

In order to describe the adiabatic quantum computer 
embedded in an environment (open system), we use the 
following rather general ansatz for the Hamiltonian 

H{t)=H,ys{t)+H,n, + Hi,-,t, (4) 

with Hsys{t) denoting the explicitely time-dependent 
Hamiltonian of the quantum computer (system) in 
Eq. |5J|, iJcnv describing the (assumed to be stationary) 
dynamics of the environment with [Hsys{t) ^ Hcnv] = 0, 
and iJint the weak interaction between them. Switching 
to the interaction picture, the total density matrix p (in- 
corporating system plus environment) evolves according 
to dtp{t) — ~i [iJint(0: (using units with h= I). For 
weak interactions between the system and the environ- 
ment governed by a small coupling constant A <C 1, it is 
quite natural to assume the following expansion of _ffi„t 

n n 
ifint = A ^ ^„ • A, + ^ <T„ • • <Tb + 0(A3) (5) 

a— 1 a, 6—1 

into one-qubit and two-qubits errors etc. The n qubits 
are labeled by a and (Ta{t) = [o'a (*)i o'a(i)i o'a(*)] the 
vector of their Pauli matrices in interaction picture, with 
the corresponding bath operators Aa{t) and Bab{t) etc. 
In the following, we shall use perturbation theory in A 
and keep only the lowest-order terms. As we shall see 
below in Eq. this amounts to omitting the Bah-term. 
(Note, however, that this approximation does not imply 
the neglect of all correlated errors.) 

Initially the system of qubits has to be realized in (or 
close to) the ground state |V'o)- Thus, we assume that 
the initial full density operator is a direct product 

p(0) = Psys(O) «) Pc„v(0) = |V'0)(V'0| «) Pcnv(O) , (6) 



i.e., system and environment are not entangled at the be- 
ginning. This should be a good approximation provided 
that the temperature is much smaller than the initial en- 
ergy gap and that the coupling i/jnt to the reservoir is 
weak. When coupling an adiabatic quantum computer to 
an environment, an essential quantity is the probability 
of measuring the ground state at the end of the com- 
putation, at time T, i.e., whether or not the algorithm 
succeeded. In the following, we shall assume perfect adi- 
abatic evolution {i.e., for A = 0, the success probability 
would be one) and hence only consider perturbations due 
to the interaction Hint , but not due to finite running time 
(nor combinations of both) . 

It turns out that the final probability of measuring the 
first excited state (Pi(T)) is a necessary and sufficient 
indicator for the success probability: The spectrum of 
the Hamiltonian JSJ consists of the ground state \Eo{t)) 
and the first excited state \Ei{t)), which come very close 
(A£;„iin = 1/Vn) at s = 1/2, whereas all iV - 2 other 
states \Ek>i(t)) have the same constant energy and are 
therefore well separated from the ground state by an en- 
ergy gap of order one. Since the temperature and hence 
the energies available in the environment are supposed 
to be much smaller than that gap of order one, tran- 
sitions from the ground state to these states \Ek>i{t)) 
are strongly (exponentially) suppressed (Pfc>i(T)) <^ 1, 
where the operators Pfe(T) = \Ek(T)){Ek(T)\ are projec- 
tors onto the fc*'' instantaneous energy eigenstate at the 
final running time T. The above observation that merely 
the ground state and the first excited state in the vicin- 
ity of the avoided level crossing (see below) are relevant 
supports our expectation that the results derived in this 
special example Q represent characteristic features of 
general avoided level crossings. 

The success of the algorithm then corresponds to 
(^0(7")) ~ 1 or alternatively (Pi(r)) < 1. For a per- 
fectly adiabatic evolution, the contributions to (Pi(T)) 
of zeroth and first order in A vanish because of Pi (T) « 
\'ipo){ipo \ and /Osys(O) = |V'o)(V'o|- For the same reason, 
the Safc-term in Eq. (|SJ do not contribute to lowest order 
©(A^). The second-order term in A gives 

(Pl(r)) ^X' J2 I'^^^l {A^a{tl)AUt2))^ 

a,h=l Q p 
fj,.u—x,y,z 

X {w^\a^{h)\w){w\<j^,it2)\w^) (7) 

after symmetrizing the integrals and neglecting terms of 
order 1 / \/N <^ 1 . This expression combines the system 
dynamics with the environment's properties: The envi- 
ronment correlation function is calculated with Ponv(O) 
and the bath operators A'^{t) and the one-qubit opera- 
tors have to be transferred to the interaction picture 

a^ait)^UUt)<Usysit), (8) 
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where Usys(t) is the unitary time evolution operator that 
is impHed by the system's Hamiltonian Hsys(t) and can 
be calculated with the known adiabatic expansion. For 
this calculation, we used the leading term only, i.e., per- 
fect adiabatic evolution and in the large- A'^ limit. 

Of the system matrix elements {w'^\a^{t)\w) only 
those with /i = x, z contribute with order one, the n = y 
term is suppressed by a factor of So, for large N 

the relevant matrix elements are 



(9) 



and the same for ('w-^\(T^(t)\w) apart from an additional 
sign (—1)™°+^, where Wa is the a"^ bit of w, i.e., the 
marked state \w) is an eigenstate of the operators with 
eigenvalues (—1)™". As one might expect, the matrix 
elements are strongly peaked at s = 1/2, i.e., at the 
point of the avoided level crossing. 

For the estimation of the time integrals in Eq. ((Jj), we 
need some information about the environment specifying 
the correlation fmiction {A'-^{ti)A'^{t2)) . In the Markov 
approximation, most of this information about the tem- 
poral behavior of the reservoir is lost. In order to com- 
pare our results with those obtained in the Markov ap- 
proximation, let us suppose that successive interactions 
are not correlated (Markovian environment) and assume 
a correlation function which is local in time 
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Insertion of this ansatz into Eq. |(7J) yields a failure prob- 
ability of 



(Pi(T))cxA2 



dt 



1 



NAE{t)J 
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for any of the three cases for s{t) discussed after Eq. 
In this case, keeping the transition errors under control 
(Pi (T)) ^ 1 requires the coupling A to be exponentially 
small with the number of qubits A — 0{1/\/N), i.e., 
the system of qubits has to be isolated exponentially 
good. Experimentally, this is not feasible and would 
therefore render the adiabatic quantum computer not 
scalable within the presence of such a Markovian envi- 
ronment. This result can be understood in the following 
way: The Markov approximation is based on the assump- 
tion that the correlations between system and reservoir 
due to their interaction at a certain point of time spread 
very fast over the entire bath and hence can be neglected 
for all interactions occurring later on. However, these 
short time-scales require excitations of comparably high 
energies. Since the adiabatic quantum computer is in 
its ground state, these high-energy excitations must have 
their origin in the bath. On the other hand, a system 
coupled to such a reservoir containing high-energy exci- 
tations for a relatively long time (the run-time of adia- 
batic quantum computation strongly increases with the 



number of qubits) is very liked to get excited - i.e., the 
computation fails. 

Instead of a Markovian bath we only assume a sta- 
tionary ([^fenvj Pcnv] = 0) rcscrvoir penv which allows for 
a Fourier decomposition of the environment correlation 
function 



(12) 



where f^l^{i^) depends on the spectral distribution (e.g., 
effective dimensionality) of the bath modes and the tem- 
perature etc. This gives 



a.b=l 



{uj)x 



g«kt+/o' drAE{T)] 1 -KO 



(13) 



plus similar terms including f^§, f^^, and with the 
associated signs (—1) and (—1)™" for x and z, respec- 
tively. Note that these terms are suppressed (in compar- 
ison to f^j^) if all qubits are coupled to the same bath, 
for example, and if we assume that the solution w has 
(nearly) the same number of zeros and ones. 

In order to evaluate the time integrations, it is use- 
ful to distinguish three different domains of u: large 
positive frequencies uj 3> A_Ei„in, small frequencies u; = 
0{AEynin), and large negative frequencies uj <C — Ai?niin- 
For large frequencies \uj\ ^ AEyain, the time integral can 
calculated via the stationary-phase (or saddle-point) ap- 
proximation. The saddle points defined by a vanish- 
ing derivative of the exponent are at w -I- AE{t'^) = 0, 
which just corresponds to energy conservation. Hence 
large positive frequencies lu ^ A£',jii,i do not contribute 
at all - which is quite natural since this corresponds to 
a transfer of a large energy to 3> Aii',jii,i from the system 
(which is in its ground state) to the reservoir. However, 
since the temperature of the bath (and hence the typical 
energy of its excitations) is typically much larger than 
AE'niin, the opposite process is possible in general. For 
small frequencies uj = 0(Ai?min), the stationary-phase 
(or saddle-point) approximation cannot be applied and 
energy conservation is also not a well-defined concept 
anymore since the rate of the external time-dependence 
is just given by 0(Aii^,„iii) and hence of the same order 
as the energies under consideration. In this case, we es- 
timate the time-integral in H13I) by omitting all phases 
(upper bound). Altogether, we get 

/+A£;„i„ 
du! f{uj) + 



2N 



duo 



AE„ 



uj^s{t*J ' 



where /(w) is understood as the appropriate sum of the 
fab^ fab^ flb^ and contributions. 
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In contrast to the first contribution, the second term 
depends on the dynamics s{t). For the three scenarios 
(s = l/T, s cx AE^, or s cx AE/y/N) discussed af- 
ter Eq. ©, the second integrand scales as Nf{—Lo)/u}^, 
f{—uj)/u!^, or y/N f{—uj)/u}^, respectively. In all three 
cases, we see that the bath modes with large frequencies 
|w| » AE^inin do not cause problems in the large- iV limit: 
The spectral function /(— w) is supposed to decrease for 
large as the bath does not contain excitations with 
large energies by assumption (low temperatures). There- 
fore, the potentially dangerous contributions stem from 
the low-energy modes w — 0{AE^in) of the reservoir, 
for which the Markov approximation is not applicable 
in general. Independent of the dynamics s{t), both, the 
first integral and the lower limit of the second integral in 
Eq. lO yield the same order of magnitude 

(P,(T)> ^ A^ffl^^. (15) 

Since AE^ain decreases as 1/Vn in the large- iV hmit, 
the spectral function /(w) must vanish in the infra- 
red limit as uj or even faster in order to keep the er- 
ror (Pi(r)) under control. Exactly at the threshold 
/(a>) ^ UJ, sub-leading contributions which scale poly- 
nomially (instead of exponentially) with the number of 
qubits n cx In Ai5,„i„ may become important: Depending 
on the reservoir (e.g., all Aa are the same versus all Aa 
are independent) , we get factors of or n from the sum 
over all qubits and for s — l/T and /(w) ^ uj, the lower 
limit of the integration yields In Ai^niin. 

Let us exemplify tha above results by means of a sim- 
ple model for the bath and consider the n qubits to be 
encoded in a chain of 1/2-spins fixed at the postions 
Ta, which are coupled to the magnetic field fiuctuations 
B{ra) via their magnetic moment 

n 
a=l 

As the reservoir, we choose a thermal photon bath (whose 
temperature is much smaller than the initial gap, but 
much bigger that the minimum gap) in D = 1,2,3 
spatial dimensions. In this case, the spectral function 
can easily be evaluated to be f[uj) ~ uj^~^ by the usual 
normal-mode expansion of the electromagnetic field into 
creation and annihilation operators. Note that, at zero 
temperature, one would gain a factor of w since the ther- 
mal particle content (Bose-Einstein distribution) has a 
1 /lv singularity in the infra-red (for a bosonic bath) . Con- 
sequently, a three-dimensional thermal photon bath does 
not destroy the scalability of the adiabatic quantum al- 
gorithm under consideration - whereas such a reservoir 
restricted to one spatial dimension (of infinite length) 
causes problems. For D = 2, the sub-leading terms scal- 
ing with poly(n) become important. 



An entirely different situation arises when the qubits 
are coupled to the amplitude of the field (representing the 
reservoir modes) instead of the field momentum density. 
This might be the case for an effectively ZJ-dimensional 
phonon bath and would result in /(w) ~ uj^~^. In such a 
situation, it is probably necessary to restrict the reservoir 
to a sufficiently finite volume (which generates an infra- 
red cut-off) in order to keep errors small. Note, however, 
that a finite volume would not improve the infra-red be- 
havior of a spin bath consisting of localized degrees of 
freedom (e.g., nuclear spins). 

In summary, we examined the effect of decoherence 
on an adiabatic quantum computer by means of Grover's 
search algorithm weakly coupled to a rather general bath 
and found that (under reasonable conditions) the infra- 
red behavior of the reservoir limits the scalability of 
the algorithm. Consequently, the Markov approximation 
does not capture the relevant features in this situation - 
instead, the spectral function /(w) of the bath provides 
a criterion to favor/disfavor certain physical implemen- 
tations. Of course, a full anal ysis would have to include 
higher orders in A (see, e.g., [l^)- If the reservoir dis- 
tributes the occupations equally between the two lowest 
levels (and there is only one level crossing), the algorithm 
still succeeds in half the cases. On the other hand, if the 
bath modes get so strongly entangled with the state |^o) 
{i.e., "stick" to it) that they prevent it from flipping to 
\w), the algorithm fails. 
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